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Abstract 



We have earlier shown that cosmic strings moving through the plasma 
at the time of a first order quark-hadron transition in the early universe can 
generate large scale baryon inhomogeneities. In this paper, we calculate de- 
tailed structure of these inhomogeneities at the quark-hadron transition. Our 
calculations show that the inhomogeneities generated by cosmic string wakes 
can strongly affect nucleosynthesis calculations. A comparison with observa- 
tional data suggests that such baryon inhomogeneities should not have existed 
at the nucleosynthesis epoch. If this disagreement holds with more accurate 
observations, then it will lead to the conclusions that cosmic string forma- 
tion scales above 10^"^ - 10^^ GeV may not be consistent with nucleosynthesis 
and CMBR observations. Alternatively, some other input in our calculation 
should be constrained, for example, if the average string velocity remains suf- 
ficiently small so that significant density perturbations are never produced at 
the QCD scale, or if strings move ultra-relativistically so that string wakes 
are very thin, trapping negligible amount of baryons. Finally, if quark-hadron 
transition is not of first order then our calculations do not apply. 
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I. INTRODUCTION 



Recent measurements of the cosmic microwave background radiation (CMBR) anisotropy 
have reached a sufficiently high level of precision that stringent bounds can be put on various 
cosmological parameters such as baryon to entropy ratio r). It is certainly quite remarkable 
that the calculations of standard big-bang nucleosynthesis (SBBN) are reasonably consis- 
tent with these stringent bounds on r/. Though several modifications to SBBN are still being 
considered to better account for the abundances of light elements. One such possibility dis- 
cussed extensively in the hterature is the so called inhomogeneous big bang nucleosynthesis 
(IBBN) [1,2] where nucleosynthesis takes place in the presence of baryon number inhomo- 
geneities. Various models have been worked out where inhomogeneitics of a particular shape 
and size are taken and their effects on nucleosynthesis are calculated. Different values of 
the light elemental abundances are calculated and compared with the observed values [2]. 
These calculations are supplemented with the investigations of baryon inhomogeneity gen- 
eration during a first order quark-hadron phase transition [3,4]. Though, it is fair to say 
that current observations do not support any strong deviation from the SBBN calculations. 
Calculations of IBBN, such as those in ref. [1,2], therefore, can be used to constrain the 
presence of baryon fluctuations in the early universe. 

In a previous paper [5], we had shown that baryon inhomogeneities on large-scales will 
be generated by cosmic string wakes during the quark-hadron transition. This arises due 
to the fact that when there are density fluctuations present in the universe, (for example, 
those which eventually lead to the formation of structure we see today), then resulting 
temperature fluctuations, even if they are of small magnitude, can affect the dynamics of a 
first order phase transition in crucial ways. There have been many discussions of the effects 
of inhomogeneities on the dynamics of a first order quark-hadron transition in the universe 
[6,7]. For example, Christiansen and Madsen have discussed [6] heterogeneous nucleation of 
hadronic bubbles due to presence of impurities. Hadronic bubbles are expected to nucleate at 
these impurities with enhanced rates. Recently, Ignatius and Schwarz have proposed [7] that 
the presence of density fluctuations (those arising from inflation) at quark-hadron transition 
will lead to splitting of the region in hot and cold regions with cold regions converting 
to hadronic phase first. Baryons will then get trapped in the (initially) hotter regions. 
Estimates of sizes and separations of such density fluctuations were made in ref. [7] using 
COBE measurements of the temperature fluctTiations in CMBR. In ref. [5], we considered the 
effect of cosmic string induced density fluctuations on quark-hadron transition and showed 
that it can lead to formation of extended planar regions of baryon inhomogeneity. 

We mention here again, as discussed in ref. [5], that there has been extensive study of 
density fluctuations generated by cosmic strings from the point of view of structure formation 
[8], and it is reasonably clear that recent measurements of temperature anisotropics in the 
microwave background by BOOMERANG, and MAXIMA experiments [9] at angular scales 
of £ ~ 200 disfavor models of structure formation based exclusively on cosmic strings [10,11]. 
However, even with present models of cosmic string network evolution, it is not ruled out 
that cosmic strings may contribute to some part in the structure formation in the universe. 
Further, cosmic strings generically arise in many Grand Unifled Theory (GUT) models. If 
the GUT scale is somewhat lower than 10^^ GeV then the resulting cosmic strings will not 
be relevant for structure formation (for a discussion of these issues, see [11]). However, they 
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may still affect various stages of the evolution of the universe in important ways. Our study 
in ref. [5] (see, also ref. [12]), and the present study are motivated from this point of view. 

In this paper we determine the detailed structure of the baryon inhomogeneities created 
by the cosmic string wakes [5] . We find that the magnitude and length scale of these inho- 
mogeneities are such that they should survive until the stage of nucleosynthesis, affecting the 
calculations of abundances of light elements. A comparison with observational data suggests 
that such baryon inhomogeneities should not have existed at the nucleosynthesis epoch. If 
this disagreement holds with more accurate observations then it will lead to the conclusion 
that cosmic string formation scales above 10^^ — 10^^ GeV may not be consistent with nu- 
cleosynthesis and CMBR observations. Alternatively, some other input in our calculation 
should be constrained, for example, if the average string velocity remains sufficiently small 
so that significant density perturbations are never produced at the QCD scale, or if strings 
move ultra-relativistically so that string wakes are very thin, trapping negligible amount 
of baryons. Of course entire discussion of this paper is applicable only when quark-hadron 
transition is of first order. 

The paper is organized in the following manner. In section II, we briefiy discuss the nature 
of density fluctuations as expected from cosmic strings moving through a relativistic fluid. In 
section III we discuss the dynamics of quark-hadron transition in the presence of string wakes, 
and discuss how baryons are concentrated in sheet like regions inside these wakes. Section 
IV discusses the results of our calculations where we present the detailed structure of the 
baryon inhomogeneities. In section V we discuss the effects, these baryon inhomogeneities 
surviving until the nucleosynthesis stage, can have on the abundances of light elements, 
and discuss the constraints on the cosmic string models arising from observations of various 
abundances. Conclusions are presented in section VI. 

II. DENSITY FLUCTUATIONS ARISING DUE TO STRAIGHT COSMIC 

STRINGS 

In this section we briefly review the structure of density fluctuations produced by a 
cosmic string moving through a relativistic fluid. The space-time around a straight cosmic 
string (along the z axis) is given by the following metric [13], 

ds^ = dt^ - dz" - dr^ - (1 - 4G'/x)V2#', (1) 

where // is the string tension. This metric describes a conical space, with a deficit angle 
of SttG/x. This metric can be put in the form of the Minkowski metric by defining angle 
ip' = {1 — 4Gfi)ilj. However, now ip' varies between and (1 — 4Gfi)27r, that is, a wedge of 
opening angle SnGfi is removed from the Minkowski space, with the two boundaries of the 
wedge being identified. It is well known that in this space-time, two geodesies going along the 
opposite sides of the string, bend towards each other [14]. This results in binary images of 
distant objects, and can lead to planar density fiuctuations. These wakes arise as the string 
moves through the background medium, giving rise to velocity impulse for the particles 
in the direction of the surface swept by the moving string. For collisionless particles the 
resulting velocity impulse is [8,15], Vimpis — ^T^GiiVst'^st (where Vst is the transverse velocity 
of the string) . This leads to a wedge like region of overdensity, with the wedge angle being 
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of order of the deficit angle, i.e. SttG// (~ IQ-^ for 10^^ GeV GUT strings). The density 
fluctuation in the wake is of order one. 

The structure of this wake is easy to sec for collisionless particles (whether non- 
relativistic, or relativistic) . Each particle trajectory passing by the string bends by an 
angle of order AnG/jL towards the string. In the string rest frame, take the string to be at the 
origin, aligned along the z axis, such that the particles are moving along the —x axis. Then 
it is easy to see that particles coming from positive x axis in the upper/lower half plane will 
all be above/below the line making an angle ^AnGfi from the negative x axis. This implies 
that the particles will overlap in the wedge of angle SnGii behind the string leading to a 
wake with density twice of the background density. One thus expects a wake with half angle 
9w and an overdensity Sp/p where [8,15], 

9^ ~ AnGfjL, ^ ~ 1. (2) 
P 

However, the case of relevance for us is cosmic strings moving through a relativistic 
plasma of elementary particles at temperatures of order 100 MeV. At that stage, it is not 
proper to take the matter as consisting of collisionless particles. A suitable description of 
matter at that stage is in terms of a relativistic fluid which wc will take to be an ideal fluid 
consisting of elementary particles. Generation of density fluctuations due to a cosmic string 
moving through a relativistic fluid has been analyzed in the literature [16-18]. The study 
in ref. [16] focused on the properties of shock formed due to supersonic motion of the string 
through the fluid. In the weak shock approximation, one finds a wake of overdensity behind 
the string. In this treatment one can not get very strong shocks with large overdensities. In 
refs. [17,18], a general relativistic treatment of the shock was given which is also applicable 
for ultra-relativistic string velocities. The treatment in ref. [18] is more complete in the 
sense that the equations of motion of a relativistic fluid are solved in the string space-time 
(Eq.(l)), and both subsonic and supersonic flows are analyzed. One finds that for supersonic 
flow, a shock develops behind the string, jiist as in the study of ref. [16,17]. In the treatment 
of ref. [18] one recovers the usual wake structure of overdensity (with the wake angle being of 
order Gp) as the string approaches ultra-relativistic velocities. Also the overdensity becomes 
of order one in this regime. 

However, it is not expected that the string will move with ultra-relativistic velocities in 
the early universe. Various simulations have shown [19] that rms velocity of string segments 
is about 0.6 for which the shock will be weak. For this case, the half angle of the wedge 
9w will also be large. We use expressions from ref. [18] which are also valid for the ultra- 
relativistic case. Resulting density fluctuation in the wake of the moving string is expressed 
in terms of fluid and sound four velocities. 



6p 167iGpuVl + u^,) 



sin9yj — — , (3) 



where Uf = Vf/^Jl — vj and Us = Vg/ ^Jl — Vg, with Vg being the sound speed. In this case, 
when string velocity Vf is ultra-relativistic, then one can get strong overdensities (of order 
1) and the angle of the wake approaches the deficit angle ~ SnGp. As mentioned in ref. [5], 
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in view of simulation results, we will use a sample value corresponding to string velocity of 
0.9 for which we take. 



0. 



w 




Sp/p ~ 3 X 10 
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(4) 



These values correspond to those obtained from Eq.3 for Vs = l/V^- In the next sec- 
tion, we will study the effects of such wakes on the dynamics of a first order quark-hadron 
transition. 

III. EFFECT OF STRING WAKES ON QUARK-HADRON TRANSITION 

In the conventional picture of the quark-hadron transition, the transition proceeds as 
follows [20]. As the universe cools below the critical temperature of the transition, 
hadronic bubbles of size larger than a certain critical size can nucleate in the QGP back- 
ground. These bubbles will then grow, coalesce, and eventually convert the QGP phase to 
the hadronic phase [20-24,6,7]. Very close to Tc the critical size of the bubbles is too large, 
and their nuclcation rate too small, to be relevant for the transition. Universe must super- 
cool down to a temperature T,c when the nucleation rate becomes significant. The actual 
duration of supercooling depends on various parameters such as the values of surface tension 
a and the latent heat L. We take the values of these parameters as in ref. [7] (motivated 
by lattice simulations [25]), a ~ 0.015T|? and L ~ 3T^. With these values, one can estimate 
the amount of supercoohng to be [7,23] (we take — 150 MeV), 



(We mention here that it has been argued in the literature that the amount of super- 
cooling may be smaller by many orders of magnitude [26]. In that situation, the density 
fluctuations of magnitude given in Eq.(4) will have even more prominent effect on the dy- 
namics of quark-hadron transition.) 

As the universe cools below Tg^, bubbles keep getting nucleated and keep expanding. This 
nucleation process is very rapid and lasts only for a temperature range of AT„ ~ 10~^, for a 
time duration of order At„ ~ lO'^tn {tn is the Hubble time) [24,7]. The latent heat released 
in the process of bubble expansion re- heats the universe. Eventually, the universe is reheated 
enough so that no further nucleation can take place. Further conversion of the QGP phase 
to the hadronic phase happens only by the expansion of bubbles which have been already 
nucleated. Even this expansion is controlled by how the latent heat is dissipated away from 
the bubble walls. Essentially, the universe cools little bit, allowing bubbles to expand and 
release more latent heat. After the phase of rapid bubble nucleation, the universe enters 
into the slow combustion phase [20]. As was discussed in ref. [5], the picture of this slow 
combustion phase is very different when cosmic string wakes are present. In the following 
we briefly review this discussion from ref. [5]. 

We take the parameters of the density fluctuations as given in Eq.(4). Density fluctuation 
Sp/p ~ 3 X 10~^ translates to a temperature fluctuation of order ATyjake = ST/T ~ 10~^. 
We note that this temperature fluctuation is larger than AT^ ~ 10~^. This means that the 
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nucleation of hadronic bubbles will get completed in the QGP region outside the wake of the 
string, while nucleation in the wake region during that stage will be suppressed. Since AT^^ 
is not too different from AT^^fee, one can conclude that the outside region will enter a slow 
combustion phase before any significant bubble nucleation can take place in the region of 
overdensity in the wake. For this it is required that the overdensity in the wake should not 
decrease in the time scale of Atn — 10~Hh- The typical average thickness of the overdense 
region in the shock dshk will be (for a wake extending across the horizon), 

dshk — sin9yj th ^ 3A;m. (6) 

Where, th = '^tn — 10 km is the horizon size, tn — 30/xsec is the age of the universe 
(at T — Tc — 150 MeV). Typical time scale, tshk: for the evolution of the overdensity in this 
region (which is smaller than the horizon size [27]), will be governed by the sound velocity 
Vg which becomes very small close to the transition temperature (e.g. ,1)5 ~ 0.1) [7,21]. We 
get, 

tsHk ^^>tH. (7) 

Vs 

Thus the density (and hence temperature) evolution in the shock region happens in a 
time scale which is too large compared to At„. It is also much larger than Attmsni— 14/isec) 
during which the quark phase is completely converted to the hadronic phase in the region 
outside the wake. We mention that in our picture, we consider the time when the universe 
has just started going through the quark-hadron transition, and we focus on a region in 
which a wake of density fluctuation of size of order rn has been created by the moving 
string. Essentially the region of study for us is the horizon volume from which the string 
is just exiting at the time when the universe temperature T = T^c- The formation of most 
of the region of shock thus happens when the temperature is still large enough compared 
to Tc so that the speed of sound is close to the value 1 / However, some portion of the 
wake will certainly form when the temperature is close enough to Tc, that the relevant sound 
speed is small, say, Vg — 0.1. The extent of the overdense region will be governed by the 
time scale tshk- Thus, if tshk is much smaller than tn, then wake will not extend across the 
horizon. 

The precise time duration, Atiag, by which the process of bubble nucleation in the shock 
region lags behind that in the region outside the wake is given by [7,24], 

Atia, ^I^fhL ^ IQ-H^^ (8) 

with Vs = l/V^- Atiag will be much larger if we take Vg = 0.1. Atiag is the extra time 
in which the temperature in the wake decreases to Tgc compared to the time when the 
temperature drops to T = Tgc in the region outside the wake. Since Atiag is at least of 
same order as Atn, we conclude that the region outside the wake enters the slow combustion 
phase before any significant bubble nucleation can take place in the wake region. It is then 
reasonable to conclude that the latent heat released in the region outside the wake will 
suppress any bubble nucleation in the wake especially near the boundaries of the wake. If 
the heat propagates to the interior of the wake, then the bubble suppression may extend 
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to the interior of the wake also, implying that there will simply be no bubbles in the entire 
wake region. In that case, hadronic bubbles which have been nucleated outside the wake 
will all coalesce and convert the entire outside region to the hadronic phase (with occasional 
QGP locahzed regions embedded in it). This hadronic region will be separated from the 
QGP region inside the wake by the interfaces at the boundaries of the wake, as shown in 
Fig.l. 

Further completion of the phase transition will happen when these interfaces move inward 
from the wake boundaries. These moving, macroscopic, interfaces may trap most of the 
baryon number in the entire region of the wake (and some neighborhood) towards the inner 
region of the wake. Finally the interfaces will merge, completing the phase transition, and 
leading to a sheet of very large baryon number density, extending across the horizon. Actual 
value of baryon density in these sheets will depend on what fraction of baryon number is 
trapped in the QGP phase by moving interfaces, and we will determine that in the following. 

It is also possible that the bubble nucleation is not entirely suppressed near the center of 
the shock region as the latent heat released by moving interfaces will have dominant effect 
locally. While the region outside the wake converts to the hadronic phase, same may happen 
at the center of the wake as well. In that case the hadronic phase will spread from inside 
the wake at the same time when the hadronic phase is moving in from outside the wake 
through the wake boundary. These two sets of interfaces will then lead to concentration of 
baryon number in two different sheet like regions, with the separation between the two sheets 
being of the order of a km or so. However, even in such a situation, the magnitude of the 
amplitude and the length scale of baryon inhomogeneity, as we determine below, will change 
only by a factor of about 2. Therefore, for simplicity of presentation, we will not consider 
this situation, and will only consider the case when only one sheet of baryon inhomogeneity 
is formed within the string wake. We also mention here that we are not considering the 
effect of density fluctuations produced by string loops. These will also lead to baryon 
number inhomogeneities via the effects discussed here. However, these structures will be on 
a more localized scale. It is more complicated to calculate the effects of density fluctuations 
by oscillating loops (especially when time scales are of crucial importance). Still, a more 
complete investigation of the effects of cosmic strings on quark-hadron transition should 
include this contribution also. 

We now determine the detailed proflle of the baryon inhomogeneity resulting from the 
above picture. For this, we have to calculate the evolution of baryon densities in the QGP 
phase and in the hadron phase as the transition proceeds. First, we note that typical 
separation between string wakes (and hence resulting baryonic sheets) will be governed by 
the number of long strings in a given horizon, which is expected to be about 15 (from 
numerical simulations [28]). The exact structure of these wakes in a given horizon volume 
needs to be known in order to study the concentration of baryons by advancing interfaces as 
the transition to hadronic phase proceeds. For example, if the string wakes are reasonably 
parallel, then they will span most of the horizon volume, as the average thickness of a wake 
(with parameters in Eq.(4)) will be order 1-2 km (a single string wake will not be expected 
to extend over the entire horizon). In such a situation, the hadronic phase will first appear 
in the regions between the wakes, which may cover a very small fraction of the horizon 
volume initially. The initial value of the fraction / of the QGP phase to the hadronic 
phase will then be close to 1. / will then slowly decrease to zero as the planar interfaces 
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(formed by the coalescence of bubbles in the region in between the overdense wakes) move 
inward, converting the QGP region inside the wake into the hadronic phase. Certainly, the 
actual situation will be more complicated than this, with string wakes extending in random 
directions, and often even overlapping. In such a situation, even the initial value of /, when 
bubble coalescence (in the regions between the wakes) forming planar interfaces, may be 
smaller than 1 (though not much smaller). However, again, this does not affect the order 
of magnitude estimate of the profile of the resulting baryon inhomogeneity as we determine 
below. Therefore, we use a simple picture, by focusing on the region relevant for only one 
string, covering about 1/15 of the horizon volume. Further, we take the initial value of / to 
be almost 1. 

The baryon evolution in this overdense wake region and outside of this region will depend 
on the detailed dynamics of the phase boundary and the expansion of the universe during 
this epoch. To study this we will follow the calculations of Fuller et.al. [22] who have 
studied the evolution of the baryon fluctuations which might have been produced at the end 
of nucleation epoch during the QCD phase transition. In ref . [22] , the evolution of baryon 
density in the QGP phase and in the hadron phase has been calculated as the hadronic region 
expands at the cost of the volume of the QGP phase during the co-existence temperature 
epoch. The main difference between their model and our model is that the QGP regions 
of interest for them are expected to be spherical, while in our case it is a thick sheet-like 
region, with planar interfaces separating the QGP region from the hadronic region. 

Let us first recall the effect of the expansion of the Universe on the dynamics of the phase 
transition [22,5]. If R{t) is the scale factor of Robertson- Walker metric, then Einstein's 
equations give [22,5], 



R{t) I STrGpit) 

where p is the average energy density of the mixed phase. The energy density, entropy 
density, and pressure {pq,Sq,Pg) in the QGP phase and {ph,Sh,Ph) in the hadronic phase 
are, 

Pq = g^aT^ + B, ^^,a^^ = ^aT^ - B (10) 

Ph = 9haT^, Sh = ^ghaT^, Ph = yaT^. (11) 

Here ~ 51 and gh — 17 are the degrees of freedom relevant for the two phases 
respectively (taking two massless quark flavors in the QGP phase, and counting other light 
particles) [22] and a = At the transition temperature we have Pq = ph which relates 
Tc and the bag constant B as, B = ^aT^{gq — gh). We define x — gqjgh to be the ratio 
of degrees of freedom between the two phases. For the mixed phase, we write the average 
value of energy density as, p = fpq + (1 — f)ph- Here / denotes the fraction of the volume 
in the QGP phase. With this, Eq.(9) can be written as. 
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Now, conservation of the energy-momentum tensor gives, 

During the transition, T and p are approximately constant. With this, Eq.(13) can be 
rewritten as 

m ^ fix -I) .^4. 

R{t) 3f{x-l) + 3' ^ ' 

Eq.(12) and Eq.(14) along with the transport rate equations, which wiU be discussed 
bellow, will give the evolution of baryon densities in the quark gluon plasma phase and in 
the hadron phase. The evolution of baryon density can be studied in each phase as the 
interfaces move forward and baryons are transported from one phase to another. If and 

are the net baryon densities in the QGP phase and the hadron phase respectively, then 
their evolution equations are as follows [22] 

nl = -nl\ + nlK-nl[^^+L-^, (15) 



= [Y3y]h<A, + nl\ + n'l^L] (16) 



where dot denotes the rate of change of the baryon density with time and Ay, arc charac- 
teristic baryon transfer rates [22] from the QGP to hadron phase and hadron to QGP phase 
respectively. The definitions of these two quantities arc discussed below. V{t) is the volume 
of the region under consideration. The term arises due to expansion of the universe and 
is given by 

m_^m (17) 

Now in our model, each cosmic string forms wake like overdensity leading the trapping of 
the QGP region in between two planar interfaces. Collapse of these two interfaces towards 
each other leads to the concentration of baryons which is the subject of study here. Numerical 
simulations have shown that in the scaling regime, there are (10 — 15) strings [28] per horizon. 
For any reasonable GUT scale, strings are definitely in the scaling regime by the stage of the 
quark hadron phase transition epoch. Initial time relevant for us is the stage when planar 
interfaces have formed (by coalescence of hadronic bubbles in the regions in between the 
wakes of different strings) at the two boundaries of overdense wakes. At this stage, we take 
the initial volume relevant for each string as, 

Vo ^ (^)r|, (18) 

where r//(= 2t) is the size of the horizon at this initial time to. Note that we take the wake 
like overdense regions to be well formed at the time to- We take the simple picture that 
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baryon concentration in each such volume is determined by the collapse of two interfaces at 
the boundary of wake of a single string, without getting significantly affected by the presence 
of wakes in the region outside the relevant comoving region. As mentioned earlier, this ap- 
proximation should be o.k. for determining the order of magnitude of baryon overdensities 
etc. Thus our representative volume is, V{t) — Vo{^-y. 

Now let us define the terms n^Ag and A/^ which appear in the transport rate equations 
Eq.(15) and Eq.(16). In our model the interface of the QGP region inside the string wake 
consists of two planar sheets. (This is in contrast to the situation in ref. [22] where the 
interface was a spherical surface). The area of each interface sheet is ^4 ~ V{ty^^, assuming 
the planar interface extending over the region with volume V{t). nlXq in ref. [22] is defined 
as the total baryon number swept by the sheets in the overdense region, and pushed through 
the underdense region, divided by the total volume in the overdense region which is f{t)V{t). 
Recall that f{t) is the fraction of the volume in the QGP phase. We get, 

„n = (19) 
"jA, ^^^^^ . [IV) 

Here, F is a filter factor which we will discuss below. (^) = is the speed of the 
interfaces. The factor 2 in the right hand side arises due to the fact there is a pair of sheets 
bounding the QGP region, which are collapsing towards each other. We write Eq.(19) as. 
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<A,= ?^^W3i^. (20) 



Similarly, n^Xh is defined as [22], 



<A, = (1)(^)(^) (21) 
- (^)(!i^)y(t)(^). (22) 

Here, Ti^ is the baryon transmission probability across the phase boundary from the 
hadron phase to the QGP phase, and Vb — {ST/mY^"^ is typical thermal velocity of baryons 
in the hadron phase, m is the mass of a nucleon. In these equations, baryon transmission 
across the interface is characterized by two parameters, F (from QGP to hadronic phase), 
and E/i (from the hadronic phase to the QGP phase). Determination of these parameters 
does not depend on the geometry of the interfaces, which is the main difference between our 
model and the one discussed in ref. [22]. For the sake of completeness, we reproduce below 
some of the steps from ref. [22] for the determination of F and S/j. 

We start with the number density of the quarks as [29], 

Uq ~ O.SgaT^, (23) 

where, g = 2nfnc is the statistical weight, and Uc, Uf are the number of colors and the 
number of flavors respectively. Following the phase-space arguments the recombination rate 
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per unit area of the quark as it approaches towards the interface separating the two phases 
has been defined in ref. [22] as, 

A = (24) 

where is the net flux of quarks and Sq is the probabihty of combining three quarks at 
the front into a color singlet, which can be estimated as [22], 

Eo = (1.4 X 10-^)SJ (25) 

where Sq characterizes the probability of transmission across the phase boundary. From 
Eq.(24) and Eq.(25) we get the total baryon recombination rate across the boundary [22] as 

A ^ (3.3 X 10^')[ Y^T.Jcm-^s-^), (26) 

where Eq.(23) has been used for the net flux of quarks. If we define ^ as the ratio of the net 
number of baryons over antibaryons to the total number of baryons, then. 



rift — ni 0.61/Xb 

,to 

h 



where the net baryon number density {riij — rii) = ^^^^^T^{y)- Therefore, the net baryon 
transport rate is given by [22] A^ = A^, i.e., 

A„ ~ (2 X 10^2)v [ 1. 112 ^ -IN (28) 

The filter factor F in Eq.(19) is defined as the ratio of the net baryon number (ANh) 
recombined to the net number of baryons encountered (Nb) at the front per unit area in 
time At. With Vz being the front velocity, the expression of A'^;, is given as 

Nb{cm-'^) = {rib - n-b)vzAt. (29) 

Similarly the expression of (ANb) can be written as 

ANb{cm~^) = AqAt 

~ (2 X 1036)EJ ]i2 X [^][^^]. (30) 

So the filter factor F is given by 

f.^ = (2.3xl0-»)|^]%.-. (31) 

So far we have considered the baryon transport rate from the QGP phase to hadron 
phase. Following the similar arguments baryon transport rate for the reverse process, i.e. 
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from the hadron phase to the QGP phase, can be calculated as follows. The net flux of 
baryons directed at the wall from the hadron phase is taken as [22], 



$'-^n^, = [^]^T^m[^]e-f, (32) 

where again m and vj, are the mass and mean velocity of a nucleon. With T^h defined above 
as the probability of a baryon to pass through the phase boundary, we can write baryon 
transport rate from the hadron phase to the quark - gluon plasma phase as [22] 

A;, « (1.1 X 10^9)[ 1—-]^ X [^]e(^)S;,. (33) 

The value of E/j depends upon the detailed dynamics of the phase boundary which can 
be calculated [30] by using chromoelectric fiux tube model. Sumiyoshi et al. [30] have shown 
that depending upon temperature this value may vary from 10^^ to 10^^ at the transition 
temperature T < 200MeK. The ratio of the two quantities S/j and can be obtained from 
the detailed balance [22] across the phase boundary for a situation when there is chemical 
equilibrium between the two phases. For this case, baryon transport rate in both directions 
are same, (i.e. Aq — Ah) , and one gets, 

^ ^ (5.4 X 10')\ ]-V^). (34) 

Using the expression for the filter factor F in terms of Eg from Eq.(31) and using 
Eqs.(19),(22), we can write the equations of the baryon transport rate (Eqs. (15-16)) in 
both regions in terms of a single parameter E/^ as follows: 



nl = 2 ^ „ " (^)[-(2.3 X 10-^)( 



x|^n^ + l.,n,Vn^[i + ^], (35) 

A^)nl - {^)v,nl] + - 3f||<. (36) 

Here Rq is the initial (when the two planar interfaces at the boundaries of the string 
wake just start collapsing) scale factor, and ^ is given in Eq.(34). 

These two equations along with the Eq.(12) and Eq.(14) have to be solved simultaneously 
to get the detailed evolution of baryon density in the trapped QGP region inside the string 
wake as well as in the region outside. Baryon inhomogeneity will be produced as baryons 
are left behind in the hadronic phase as the interfaces collapse. We now study the profile 
of the resulting baryon overdensity after the interfaces collapse away. Let Nq{t) be the total 
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baryon number in the QGP region at a particular time t. Nq{t) is related to the baryon 
density nl as : 

N,{t)=nlit)Vit)fit). (37) 

Taking center of the wake as the origin and considering motion of the interface along z 
direction, we can write 

f{t)V{t)^2A{t)zit). (38) 
From this we get the evolution of the thickness z{t) as, 

. Mv.<i>^. (39) 

To get the profile of the baryon inhomogeneities, let p{z) be the baryon density which is 
left behind at position z as the interfaces collapse. We get, 

Ng{z) -Ng{z + dz) = Adzp{z), (40) 

where the time dependence of z is given in Eq.(39). We get. 

Thus we finally get the density of baryons p{z) left behind by the collapsing interfaces 

as. 

Note here that derivation of this equation assumes that baryons left behind by the 
collapsing interfaces remain in the same region of z, and do not diffuse away. On the other 
hand, the derivation of equations for baryon transport (Eq.(35)-Eq.(36)) was based on the 
assumption that baryons in both phases homogenize, so that baryon transport equations 
can be written only in terms of two baryon densities, one for each phase. If baryons do 
not homogenize in the hadronic phase (as was assumed for Eq.(42)), then it will increase 
the reverse baryon transport rate, i.e. from the hadronic phase to the QGP phase. This 
will only increase the baryon inhomogeneity produced. Also, as mentioned above, values 
of Eft are expected to be very small. We find that even with two orders of magnitude 
increase in the value of S^, the relevant width of the profile of baryon overdensity p{z) only 
increases by one order of magnitude. Thus, within this uncertainty, we will use Eq.(42) 
to determine the baryon inhomogeneity profile. Finally we mention that baryon diffusion 
length for the relevant overdensities here always remains less than few cm, while the length 
scales of inhomogeneities of interest to us are at least one order of magnitude larger. 
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IV. RESULTS 



Eqs.(12),(14),(35),(36) are numerically solved simultaneously to get the evolution of 
baryon densities n^, and in the two phases for two different values of critical temper- 
ature Tc = 150 MeV and = 170 MeV, and for two values of = 10"^ and = 10"^ 
Figs. 2-3 show plots of nl for these cases. Resulting profiles of baryon ovcrdcnsity p{z) are 
calculated using Eq.(37) and Eq.(42), and are shown in Figs. 4-5. Note that in Fig.5 there 
are wiggles in the plot of p{z). This is due to numerical errors in calculating Ng{t). As ng{t) 
increases, f{t)V{t) decreases, leading to extremely slow variation in Nq{t) for most of the 
time duration. Thus, as variations in nq{t) and in f{t)V{t) compensate for each other, the 
errors become relatively large, as seen in Fig.5. We have checked that these errors are in 
better control for other parameters where the change in Nq{t) is larger. For example, change 
in Ng will be expected to be larger if Eg (which determines baryon transport rate from the 
QGP phase through the interface to the hadronic phase) is made larger while keeping S/j 
(determining baryon transport from hadronic to QGP phase) fixed. This could be achieved 
by taking the nucleon mass m in Eq.(34) to be smaller. We have verified that indeed this 
happens. For smaller values of m errors become much lower. Of course, as m is not a free 
parameter, we do not give plots for different values of m. Also, for most of time duration, 
baryon fiow from the QGP phase to hadronic phase dominates over the reverse fiow. With 
m fixed, when E/j is increased. Eg also increases proportionally via Eq.(34). Therefore, a 
larger E^ again leads to a more rapid variation of N^, giving better control of errors. This 
can be seen from plots in Fig.4 which correspond to Hh = 10^^. No wiggles are seen here, 
compared to the situation in Fig.5 which corresponds to E/i = 10^'^. 

We have used Mathematica routines for numerically solving these coupled differential 
equations. Due to very wide range of numerical values of various parameters involved, time 
step for solving differential equation had to be chosen judiciously. For example, for initial 
times, when distance scale of region in between the interfaces is about 1 km, large values of 
time step is chosen. As the distance scale decreases, the time step is decreased by factors of 
100, ranging from 0.1 pscc to 10^^/xsec. This gives a good overall control on the accuracy. 
An indicator for the error in the numerical solution we obtain is the value of total baryon 
number A^^ = Ng + N^. As the interfaces collapse, converting the QGP phase to the hadronic 
phase, Ng decreases while Nh increases. However, Nt must remain constant. We find that 
the value of Nt remains reasonably constant over the entire range of integration relevant to 
us (as shown in Figs. 4, 5). There is a tendency of small net increase in Nt as a function of 
time. The net increase in the value of Nt (which indicates error in the numerical evolution) 
remains less than 5% of the net change in the value of Ng over the range of integration in 
the plots. From Eq.(42), we see that the value of p{z) is directly proportional to dNg/dz. 
Only time dependence in the proportionality factor is in -R(t), which changes little over the 
entire time period, and its evolution is smooth, without any random errors. Similarly, the 
evolution of z{t) is smooth, without any random errors. Thus, resulting error in p{z), as 
shown in the plots below, should also be less than about 5%. Apart from this error, there is 
also a random component in the error (again, resulting from extremely slow variation of Ng), 
leading to random wiggles in the plots of p{z) as visible in Fig.5. The largest magnitude of 
this error in p[z) is about 20%. (This error is negligible for plots in Fig.4, and also much 
smaller than 5% at other parts of plots in Fig.5 where wiggles are not seen.) As our interest 
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is in order of magnitude estimates of the baryon overdensity (its magnitude, as well as its 
spatial profile), even the largest possible error of 20% here does not affect our results and 
conclusions. 

As we will discuss in the next section, relevant values of the overdensity n'^/ni, for us 
is about 1000. Here and ni, are baryon densities in the overdense and the background 
regions respectively. From above plots wc sec that for = 10^^, the thickness of the region 
inside which n'^^/nh > 1000 is about 5 meters for Tc = 150 and about 4 meters for Tc = 
170 MeV. For E/j = 10~^ this thickness varies from about 0.5 meters to about 4 meters as 
Tc changes from 170 to 150 MeV. As baryon density sharply rises for small z, it is more 
appropriate to calculate the largest value of the width of the inhomogeneity region within 
which average value of baryon density is 1000 times larger than the asymptotic baryon 
density. We find that this width is at least an order of magnitude larger than the values 
mentioned above. For E/j = 10~^ this width is about 100 m for Tc = 150 MeV and is about 
60 m for = 170 MeV. For E^ = 10"^ the values of this width are about 120 m and 90 
m for Tc = 150 MeV and 170 MeV respectively. As we will see below this type of baryon 
overdensities can strongly affect abundances of light elements, thereby constraining various 
parameters of cosmic string models. 

V. NUCLEOSYNTHESIS CONSTRAINTS 

With the baryon inhomogeneity profile determined as above at the QCD scale, we need to 
know the amplitude and length scale of this inhomogeneity at the epoch of nucleosynthesis. 
For this we use results of ref. [4] where evolution of baryon inhomogeneities of varying 
amplitude and length scales has been analyzed in detail. From ref. [4] one can see that 
baryon inhomogeneities of initial magnitude n[,/nf, ~ 1000 at the QCD scale should survive 
relatively without any dissipation until the stage when temperature T ~ 1 MeV for all the 
values of length scales relevant for us, i.e. few tens of cm upto about 100 meters. (For 
example inhomogeneities with baryon to entropy ratio of about 10~^ almost do not change 
during their evolution. Inhomogeneities with larger amplitude eventually dissipate to this 
value. See, ref. [4].) Though, the length scales in ref. [4] are taken to be comoving at 100 
MeV, the results there should apply relatively unchanged for the the values of Tc we have 
considered, i.e. Tc = 150 and 170 MeV. 

To study the effects of these resulting inhomogeneities at the nucleosynthesis epoch, we 
use the results of the calculations of the IBBN model developed by Kainulainen et al. [2]. 
The four crucial parameters in this model are, the average baryon density {riavg)^ the density 
contrast (i? = ri'i^/rib), the volume fraction of the high density region, and the distance 
scale r of the inhomogeneity at the onset of nucleosynthesis. Out of these, the last three 
parameters characterize the properties of baryon inhomogeneity regions. We obtain these 
three parameters from our model and check them with the numerical results in ref. [2] to 
determine their effects on nucleosynthesis results. Though, the geometry in our case is not 
exactly the same as the spherically symmetric geometries that have been considered in ref. 
[2]. However, note that the planar sheet like inhomogeneities of our model should be similar 
to the geometry of spherical shell (SS) considered in [2]. Therefore, for rough estimates, we 
will simply take the results of inhomogeneities of the shape of spherical shell from ref. [2] , 
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and apply it to our case, making sure of using the corresponding values of the parameters 

R, fy, and r. 

The results in ref. [2] for the SS geometry were given for a fixed value oi R = 1000, with 
fy varying from about 0.023 to 0.578. In order to be able to use the results of ref. [2], we 
therefore determine the thickness (and hence the value of f^) of the baryon inhomogeneity 
regions from Figs.(4)-(5) within which R > 1000. Again, note that the plots in Figs.(4)-(5) 
are given for baryon inhomogeneitics at the QCD scale. However, results of ref. [4] show that 
there is no significant dissipation of these inhomogeneitics upto the nucleosynthesis scale. 
Thus, with length suitably scaled with the scale factor of the universe, profiles in Figs. (4)- 
(5) can be used at the nucleosynthesis stage. We note from Fig. 4, for E/j = 10~^ that the 
region of baryon inhomogeneity within which R > 1000 has a value of ~ 5/2000 ~ 0.0025 
for both values of Tc = 150 and 170 MeV. Here the relevant size of the whole region is 
taken to be about 2 km. This value of fy is smaller than the smallest value of fy ~ 0.023 
considered in ref. [2] for the SS geometry case. However, above estimates for fy are clearly 
an underestimate as the baryon density is sharply peaked inside the overdense region. As 
discussed earlier, if we calculate the largest value of the width of the inhomogeneity region 
inside which the average value of the baryon density is 1000 times larger than the asymptotic 
value then the resulting widths are very large, varying from about 60 meters to about 100 
meters. This will then lead to a value of fy of about 0.03 - 0.05 which are sufficiently large. 
Note also, that crucial parameter for our case, using which one can determine the order of 
magnitude effects of baryon density fluctuations on element abundances, is the optimum 
value of the parameter r. This value depends very weakly on fy for the SS geometry, with 
Topt ~ fv^^^ for /i) << 1 (see, ref. [2]). Thus, even with smaller estimates of fy, the value of 
Topt relevant for our case will be only about factor 2 larger than the value in ref. [2] for the 
case fy ~ 0.023. Similarly, from Fig.5 for the case E^j = 10~^, we see that thickness of the 
inhomogeneity within which R > 1000 is about 4 m for Tc = 150 McV, and about 0.5 m for 
Tc = 170 MeV. Corresponding values of fy arc about 10~^ and 10^^ respectively. In these 
cases, value of Vopt will be increased by about one order of magnitude. Again note that if we 
take the average baryon density then the relevant width is much larger, about 90 meters to 
120 meters. This then leads to a large value of fy, about 0.045 to 0.06, and hence estimate of 
Topt remains unchanged. Note also, that in ref. [2] it is mentioned that for maximum effect, 
the value of Rfy should be much greater than about 7 (the SBBN proton/neutron ratio at 
the onset of nucleosynthesis). The smallest value of Rfy considered in ref. [2] is 23. In our 
case for smaller estimates of fy the value of Rfy is 2.5. However, when we take average 
baryon density, then the value of Rfy ranges from about 30 to 60 which is similar to the 
values considered in ref. [2]. 

Next thing we note is that the typical separation between the inhomogeneitics, i.e. the 
parameter r, is about 1-2 km for our case. This corresponds to about 100- 200 km length 
scale at the nucleosynthesis epoch. Importantly, this is precisely the range of values of r for 
having optimum effects on nucleosynthesis calculations in ref. [2] . Even with the variations 
in fy as discussed above, one can conclude that with R = 1000, and with values of /„ 
corresponding to different cases in Figs. (4)- (5), the length scales of inhomogeneitics in our 
model (inter-inhomogcncity separation) is roughly in the right range to have optimum effects 
on nucleosynthesis calculations. 

We now apply observational constraints on the abundances of various elements. The 
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most basic constraint is on the abundance of ^He by mass, denoted by Y. If we take a 
liberal range of values oi Y = 0.238 — 0.244 (sec, rcf. [9]), then using the results of IBBN 
calculations in ref . [2] , we see that for inhomogeneities with optimum value of r (which we 
have argued to be the case here), the corresponding value of rj is between about 4 x 10~^° 
to about 8 X 10~^°. (We mention that most plots in ref. [2] have been given for a different, 
centrally condensed geometry of the inhomogeneities. However, it has been mentioned there 
that for SS geometry also results are not too different.) These values are about a factor 
2 larger than the allowed values of rj for the case of SBBN. (Abundance of ^Li for IBBN 
models favors somewhat smaller values of rj. One needs a careful and detailed comparison of 
abundances of various elements, ^He, ^Li, and D. However, in view of various uncertainties 
of our model we will only consider the case of ^He here.) 

An independent estimate of rj comes from the cosmic microwave background (CMBR) 
anisotropy measurements. Constraints coming from various experiments seem to constrain 
rj to be less than 6 x 10~^*^. If one takes large estimates of ^He, then IBBN calculations 
suggest that the corresponding value of rj will not be consistent with the value obtained 
from CMBR measurements. Note that SBBN estimates of 7] for the above range of Y are 
in very comfortable agreement with CMBR measurements. With this, we conclude that 
it is suggestive that the baryon inhomogeneities of the type produced by cosmic strings 
as discussed above are not consistent with the combined observations of '^He abundance 
and CMBR anisotropy measurements. Therefore, some of the parameters of the cosmic 
string model may have to be constrained, so that such inhomogeneities are not produced 
at the QCD scale. Of course, this is assuming that quark-hadron transition is a first order 
transition. If the transition is of second order, or a cross-over, (as suggested by many 
studies), then our calculations do not apply. 

If the transition remains of first order, then there are several ways in which production 
of such inhomogeneities can be avoided. First, if the value of string scale is smaller, say 
~ 10^^ GeV, then from Eq.(3) we see that resulting value of Sp/p will be smaller by one 
order of magnitude. This implies that the excess temperature inside the string wake region 
ST/T will be about 10~^. This value is much smaller than the value of ATg^ required for the 
supercooling for bubble nucleation to start in the outside region. In such a situation, bubble 
nucleation inside the wake will not be entirely suppressed, though it may still lag behind 
the nucleation of bubbles in the outside region. Therefore, it is still not excluded that some 
sort of large scale baryon inhomogeneities will get produced even with string scale of 10^^ 
GeV. If this string scale was smaller than 10^"^ GeV, then resulting value of 6T/T will be 
even smaller than AT„ ~ 10~^. It is extremely unlikely that in such a case any significant 
effect will be there on the dynamics of quark-hadron phase transition due to the presence 
of string wakes. 

Yet another possibility is that string velocity is cither much smaller, or extremely close 
to the speed of light. In the first situation, resulting value of Sp/p is very small, so no effect 
will be there on the transition (just as the case for small G/j). For the second situation, 
when strings move ultra- relativistically, Sp/p will be very large (of order 1), so quark-hadron 
transition dynamics will be strongly affected, producing sheet like baryon inhomogeneities. 
However, in this case the wake angle 6^ will be very small (of order SnGp). In such a 
situation, for the region relevant for a single cosmic string, string wake will cover a very 
small fraction of the total volume. Thus, when the region outside the wake undergoes 
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hadronization, many localized regions of baryon inhomogeneities will form just as in the 
conventional models of homogeneous nucleation. Planar interfaces will still form, but they 
will be able to concentrate baryons from only a very small fraction of the total volume. 
Thus, resulting baryon inhomogeneities will contribute to negligible baryon fluctuation on 
the average. 

In summary we conclude that observational constraints from abundance of ^He, and 
CMBR anisotropy measurements may constrain the cosmic string scale to be less than 
10^'' — 10^^ GeV. Though there are many uncertainties in our model. Alternatively, string 
velocities, either need to be very small, or very large. Though as strings are in not in the 
friction dominated regime (for relevant values of GUT scale), it may be harder to decrease 
average string velocity sufficiently. Similarly, due to random velocity components of different 
string segments, it may not be easy to argue for ultra-relativistic motion of strings. Of course 
entire discussion of this paper (as well as in ref. [5]) is applicable only when quark-hadron 
transition is of first order. 

VI. CONCLUSION 

We have calculated the detailed structure of the baryon inhomogeneities created by the 
cosmic string wakes [5] . We find that the magnitude and length scale of these inhomogeneities 
is such that they survive until the stage of nucleosynthesis, affecting the calculations of abun- 
dances of light elements. A comparison with observational data suggests that such baryon 
inhomogeneities should not have existed at the nucleosynthesis epoch. If this disagreement 
holds with more detailed calculations and more accurate observations, then it will lead to 
the conclusion that cosmic string formation scales above a value of about 10^^ — 10^^ GeV 
are not consistent with nucleosynthesis and CMBR observations. Alternatively, some other 
input in our calculation should be constrained, for example, the average string velocity can 
be sufficiently small so that significant density perturbations are never produced at the QCD 
scale, or strings may move ultra-relativistically so that resulting wakes are very thin, and 
trap a negligible amount of baryon number. Finally, all these considerations are valid only 
when quark-hadron transition is of first order. 

There are many uncertainties in our model, for example treatment of multiple wakes is 
rather ad hoc. Similarly, we have tried to use results from ref. [2] adopting them for our case 
even though detailed geometry of baryon inhomogeneity in our case is different. A more 
careful, detailed calculation of abundances of elements is needed for the present case. 

The uncertainties in various observations of abundances of elements, as well as CMBR 
anisotropy will be reduced as precision of various measurements gets better. Then one will 
be able to say with a greater certainty whether IBBN results puts a strong restriction on the 
density ffuctuations, and hence on cosmic string parameters, or the order of quark-hadron 
phase transition. 
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FIG. 1. 

This figure shows a portion of the overdense shock region. Q and H denote QGP and 
hadronic phases respectively and sohd hnes denote the interfaces separating the two phases. 
Shaded region, marked as B, denotes resulting sheet like baryon inhomogeneity. 
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FIG. 2. 

These figures show plots of baryon density nl in the QGP phase inside the wake, as a 
function of time. S/j here is 10~^. Top figure is for = 150 MeV and bottom figure is for 
Tc — 170 MeV. Here and in Fig. 3, nl is in fm~^, while time is given in /is. 
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FIG. 3. 

Plots of nl as a function of time. E/j here is 10~^. Top figure is for T^, — 150 MeV and the 
bottom figure is for Tc = 170 MeV. 
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FIG. 4. 

These figures show profiles of baryon inhomogeneities p{z) generated by collapsing planar 
interfaces. E/j here is 10~^. Top figure is for Tc = 150 MeV and the bottom figure is for 
Tc — 170 MeV. Here, and in Fig.5, p is in units of fm"^ while z is given in meters. Insets 
show expanded plots of the region where p becomes larger than 1000 times the asymptotic 
value. We have estimated the error in numerical evaluation of p{z) (here, and in Fig.5). 
Largest error is about 20 % and occurs whegg wiggles are seen in Fig.5. At other parts of 
plots, error remains below about 5%. 
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FIG. 5. 

Plots of p{z) vs. z for the case when S/j is 10~^. Top figure is for Tc = 150 MeV and the 
bottom figure is for = 170 MeV. 
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